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►^ I The CRESST-II collaboration have announced evidence for the direct detection of dark 

C^- ■ matter in 730 kg-days exposure of a CaW04 target. We examine these new results, along 

^ ■ with DAMA and CoGeNT data, in the context of the mirror dark matter framework. 

CSl '• We show that all three experiments can be simultaneously explained via kinetic mixing 

cri \ induced elastic scattering of a mirror metal component off target nuclei. This metal 

^ \ component can be as heavy as Fe' if the galactic rotational velocity is relatively low: 

Vrot ~ 220 km/s. This explanation is consistent with the constraints from the other 
experiments, such as CDMS/Ge, CDMS/Si and XENONIOO when modest ~ 20 - 30% 
L^ . uncertainties in energy scale are considered. 
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1 Introduction 

Over the last decade or so, progress has been made in efforts to directly detect darlc mat- 
ter. In particular, the DAMA/Nal [1] and DAMA/Libra [2] experiments have obtained 
very exciting results in their dark matter search. Recall that these experiments have 
observed a modulation in the 'single hit' event rate with a period and phase consistent 
with expectations from dark matter interactions [3]. Background rates are expected to 
be time independent with the possible exception of muon induced backgrounds. How- 
ever it has been known for a long time that muons cannot mimic the dark matter annual 
modulation signature [4]. The DAMA experiments thus provide very convincing evidence 
that dark matter interactions have been detected. 

Another interesting development is the results obtained in the CoGeNT experiment [5]. 
That experiment features a Germanium target and very low energy threshold. The Co- 
GeNT team observed a rising event rate at low energies. These events could not be 
explained by known backgrounds and can be interpreted as evidence supporting the 
DAMA signal. Very recently, the CRESST-II collaboration have announced results for 
their dark matter search with 730 kg-days of net exposure in a CaW04 target [6]. The 
CRESST-II data also cannot be explained by known backgrounds and are compatible 
with a dark matter signal also rising at low energies. 

Attempts to explain the positive signals of DAMA, CoGeNT and CRESST-II in terms 
of standard WIMP dark matter have proven challenging [7]. However, it has been known 
for some time[8, 9, 10], that mirror dark matter has a number of novel and distinctive 
features which make it a suitable candidate to explain the data. In particular mirror 
dark matter features a mass dependent velocity dispersion vq oc J^ , Er dependent 

Rutherford scattering ^^- oc -^, and particles necessarily in the relevant low mass range. 

With these features the DAMA and CoGeNT data can be simultaneously explained[ll, 
12]. The main purpose of this paper is to examine the implications of the new CRESST- 
II results for mirror dark matter. In the process we will update the CoGeNT analysis 
taking into account their recently reported surface event correction factor [13]. It turns 
out that these developments point to new parameter regions within the mirror dark 
matter framework. Heavy mirror elements, ~ Fe', can explain all three experiments if 
the galactic rotational velocity is low Vrot ~ 220 km/s. Other regions of parameter space 
are possible. We explore a few examples in some detail and examine the question of 
their compatibility with constraints from sensitive but high threshold experiments such 
as XENONIOO and CDMS/Ge. 

2 Mirror dark matter and its direct detection 

Mirror dark matter posits that the inferred dark matter in the Universe arises from a 
hidden sector which is an exact copy of the standard model sector [14] (for a review and 
more complete list of references see ref. [15] )^. That is, a spectrum of dark matter particles 
of known masses are predicted: e', H', He', O', Fe',... (with iriei = nie, mn' = itlh, 



^Note that successful big bang nucleosynthesis and successful large scale structure requires effectively 
asymmetric initial conditions in the early Universe, T' <^ T and rib' /rib ~ 5. See ref. [16] for further 
discussions. 



etc). Kinetic mixing of the f/(l)y and its mirror counterpart allows ordinary and mirror 
particles to interact with each other[17]. This f/(l)y kinetic mixing induces photon- 
mirror photon kinetic mixing: 
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where F^y is field strength tensor for the photon and F' is the field strength tensor for 
the mirror photon. This interaction enables charged mirror sector particles of charge e 
to couple to ordinary photons with electric charge ee [18]. The cross-section of such a 
particle, say a mirror nucleus, A' , with atomic number Z' and velocity v to elastically 
scatter off an ordinary nucleus (presumed at rest with mass and atomic numbers A, Z) 
is given by[8]:" 
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where 



A = Fl{qrA)Fl{qrA') 3 

niA 

and FAiqTA) [-Fa'(?^A')] is the form factor which takes into account the finite size of 
the nucleus [mirror nucleus]. A simple analytic expression for the form factor, which we 
adopt in our numerical work, is the one proposed by Helm[19, 20]. 

What about the distribution of mirror particles in the galactic halo? These are 
assumed to be, predominately,^ roughly spherically distributed in a multi-component 
plasma containing e', H', He', O', Fe',...[23]. The interaction length is typically much less 
than a parsec[24] and the dark matter particles form a pressure-supported halo, ideally 
described by a common temperature, T ^. This temperature can be estimated from the 
condition of hydrostatic equilibrium which implies: 

T = \mvl, (4) 

where Vrot is the galactic rotational velocity and ffi = J2''^a'1tia'/ J2''^a' is the mean 
mass of the particles in the halo. The halo distribution of mirror species is then J'a' = 
exp{—E/T) = exp{—^mA'V^/T) = exp{—v'^/vl). Evidently 
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■^Unless otherwise indicated, we employ natural units where Ti — c—1. 

^There may also be a remnant disk of dark mirror stars [21], but constraints from gravitational 
microlensing (MACHOs searches) limit this component to be subdominant[22]. 

^There are limits on self interactions of dark matter from observations of the Bullet cluster [25]. 
This sets stringent limits on self interactions provided that the bulk of the dark matter particles are 
distributed throughout the cluster and not bound to individual galaxies. However mirror dark matter 
is dissipative and in clusters (or at least in some of them) the bulk of the dark matter particles might 
be confined in galactic halos [c.f. ref. [26]]. Under this assumption mirror dark matter is consistent with 
Bullet cluster observations. 



The parameter m can be estimated within the mirror dark matter modeL For a H', 
He' mass dominated halo, taking into account that the plasma is expected to be fully 
ionized, we have 
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where C,He' is the mass fraction of He'. The primordial He' abundance can be computed 
as a function of the kinetic mixing parameter, e, and for e ~ 10^^ the calculations suggest 
that the primordial He' mass fraction, Y^, is around 0.9 [27]. If C,He' ~ ^' then Eq.(6) 
suggests that m ^ 1.1 GeV. Even if mirror stellar evolution depleted H' further (or 
somehow depleted both H' and He'), this would have only a relatively minor affect on 
the estimation of rh. 

A key feature of this model is that the heavy components with rriA' ^ fn have 
Vo[^'] <^ Vrot- This can drastically reduce the tail of the distribution [c.f. standard 
WIMPs which have Vq = Vrot]- This feature can help explain why higher threshold 
experiments such as CDMS/Ge and XENONIOO currently do not see a signal while the 
lower threshold DAMA and CoGeNT experiments do. 

The differential scattering rate for A' elastic scattering on a target nuclei. A, is given 
by^: 

NtUa' / -77; ; \\\d V (7) 
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where the integration limit, Vmin, is given by the kinematic relation: 



{mA + mA'YER 
\ 2mAm\, ■ ^^^ 

In Eq.(7), Nt is the number of target nuclei per kg of detector and Ha' is the number 
density of halo dark matter particles. A' , at the Earth's location. This number density 
can be expressed in terms of the halo mass fraction of species A' , ^a'-, and total mass 
density, pdm via n^' = Pdm^A' /'rriA' (we fix pdm = 0.3 GeV/cm ). Also in Eq.(7) v is the 
velocity of the halo particles relative to the Earth and v^; is the velocity of the Earth 
relative to the galactic halo^. The halo distribution function, in the reference frame of 
the Earth, is then given by the Maxwellian distribution: 

The integral, Eq.(7), can easily be evaluated in terms of error functions[10, 20] and 
numerically solved. 



^The upper limit of integration in Eq.(7) is taken as infinity since we are dealing with dark matter 
particles with significant self interactions. The self interactions will prevent particles in the high velocity 
tail of the Maxwellian distribution from escaping the galaxy. 

^In all numerical work we include an estimate of the Sun's peculiar velocity so that {ve) = Vrot + 12 
km/s. 
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10.2 - 13.0 
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3.2 


13- 16 
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6.1 


16- 19 


11 


7.0 
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11.5 


25-40 
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20.1 



Table 1: CRESST-II data: total number of events and estimated background. 

To compare with the measured event rate, we must include detector resolution effects 
and overall detection efficiency, ef{E^), when the latter is not already included in the 
experimental results: 

^^ ef{E^)^ /i^e-(^-~^H)V2<x?..rf^^ (10) 



dE"^ '' "" ^2^ ares J dE, 

where E^ is the measured energy and (Tres describes the resolution. The measured energy 
is typically in keVee units (ionization/scintillation energy). For nuclear recoils in the 
absence of any channeling, keVee = keV/g, where g < 1 is the relevant quenching factor. 
Channeled events, where scattered target atoms travel down crystal axis and planes, 
have g ~ 1. In light of recent theoretical studies [28], we assume that the channeling 
fraction is negligible. Channeling, though, is a complicated subject and the modeling 
used in ref.[28] might miss important physics. Therefore it is, of course, still possible 
that channeling could play an important role, which could modify the favored regions of 
parameter space somewhat. 

3 Mirror dark matter interpretation of CRESST-II, 
DAMA/Libra and CoGeNT signals 

The CRESST-II experiment 

The CRESST-II data[6] arises from 8 detector modules, with recoil energy thresholds 
(keV) of 10.2, 12.1, 12.3, 12.9, 15.0, 15.5, 16.2, 19.0. The data is plotted in figure 1. With 
a view to performing a standard x^ analysis of this data, we bin this data into 5 bins, 
summarized in table 1. This table also indicates the expected background rate estimated 
from all known sources of background [6]. 
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Figure 1: CRESST-II data. 



The rates, Ri, relevant for the CRESST-II experiment are computed using Eq.(7) and 
Eq.(lO). As discussed earUer, the halo distributions in Eq.(7) depend on the parameter 
m which has been estimated to be m ?a 1.1 GeV. A reasonable uncertainty in rh is ±20%, 
but such a variation does not significantly affect any of our conclusions. We therefore 
fix m = 1.1 GeV in all of our numerical work. We use ares = 0.3 keV[6] and assume 
detection efficiencies for the three target elements: e/ = 0.90 for O, W and ej = 1.0 for 
Ca, which take into account their acceptance region. As per previous analysis [11, 12] we 
make the simplifying assumption that the mirror metal components are dominated by 
a single element A' in a given experiment. This can only be an approximation, however 
it can often be a reasonable one given the relatively narrow energy range probed in the 
experiments [e.g. the signal region in the experiments are mainly: 2-4 keVee (DAMA), 
0.5-1 keVee (CoGeNT), 12-14 keV (CRESST-II)]. With this assumption and assuming 
a particular value for Vrot, the fit to the data depends on two parameters rriA', ^V^- ^ 
We define x^ for the CRESST-II data in the usual way: 



X^{mA',e^JU') = Y1 
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where Ri is the theoretically predicted rate and Bi is the estimated background in the 
jth g^ergy bin. We evaluate x^ as per Eq.(ll), with the constraint rriA' < rriFe — 55.8mp. 
We have not considered any uncertainty in the CRESST-II energy scale. 



^In our numerical work we allow A' , Z' to have non-integer values, with Z' = A' /2 (except when we 
specifically consider A' = Fe', then Z' = 26 and A' ~ 56). Since the realistic case involves a spectrum 
of elements, the effective mass number can be non-integer. 



The DAM A annual modulation signal 

For DAMA we analyse the annual modulation signal using the 12 bins of width 0.5 
keVee in the energy range: 2 — 8 keVee[2] taking into account the detector resolution [29]. 
An important issue is the nuclear recoil energy scale which is set by the values of the 
quenching factors gTVa, qi- In the recent analysis[ll, 12] we used q^a = 0.30 ± 0.06, 
qi = 0.09 ±0.02, where the central values were taken from averages over a large range of 
energies [1]. However quenching factors are in general energy dependent and many mea- 
surements of qi tend to show increasing values in the important low energy region. This is 
also supported by some theoretical arguments based on the classical Birks forumula[30]. 
[See e.g. ref. [31] for a summary figure of much of the available information]. In view of 
the current situation, we allow for a range of possible values for the quenching factors: 

g^a = 0.28 ±0.08, g/ = 0.12 ± 0.08 (12) 

We minimize x^ varying g^va, qj over the above range of values (assumed energy inde- 
pendent for simplicity). Values of qj, q^a outside the above range are certainly possible, 
such as the higher values qNa ~ 0.6, qi ~ 0.3 suggested by Tretyak[30]. While we don't 
specifically consider the possibility of such high quenching factors in the numerical work 
presented here, qualitatively we note that higher quenching factors generally move the 
DAMA allowed region to lower values of e^J^A'-, f^A' and also yields interesting parameter 
space overlapping with that from CoGeNT and CRESST-II. 

The CoGeNT experiment 

Our analysis procedure for CoGeNT is different to our earlier analysis[ll, 12]. It 
has been pointed out that there may be surface events contaminating CoGeNT's signal 
region[13] which have not been excluded by CoGeNT's rise time cut. Preliminary es- 
timates indicate that this can be a very large effect - the signal may be reduced by a 
factor of around 0.3 in the important low energy region [E^i < 1.0 keVee)[13, 32]. In 
view of this we analyse CoGeNT data using 15 bins of width 0.1 keVee in the region 
0.5 — 2.0 keVee. We take the efficiency corrected and surface event corrected CoGeNT 
data from figure 4 of ref. [32]. We allow a constant background, which we fit to the data 
in this energy range. We have taken into account some uncertainties in energy scale by 
minimizing the x^ for CoGeNT over the variation in quenching factor, qce = 0.21 ±0.04. 
Obviously, the CoGeNT spectrum is quite uncertain at the present time and one could 
easily assign ~ 40% uncertainty in the rate. Such an uncertainty would imply a ~ 20% 
uncertainty in inferred values of e^/^ from CoGeNT data. This is in addition to the 
uncertainties we have considered. Thus, the reader should be aware that the favored 
region we derive for CoGeNT might move considerably as future data is accumulated 
and the backgrounds are better understood. 

The analysis 

Table 2 gives the Xmin values obtained from the fit to the DAMA annual modulation 
and CoGeNT, CRESST-II spectrum data for a fixed value of Vrot = 200 km/s, which we 
have found as an example where all three experiments have overlapping favored regions 
of parameter space. The value Vrot = 200 km/s is around 20% less than some recent 
estimates of this quantity [33]. However allowing for a small bulk halo rotation, e.g. in a 
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Experiment 



X^{iTLin)/d.o.f. 



Best fit parameters 



DAMA (annual mod.) 5.7/10 

CoGeNT (spectrum) 15.7/12 

CRESST (spectrum) 2.6/3 



nip 


= 55.8 , 


eVU/lO-'' = 2.5 




= 42.0 , 


eVU/lO-'' = 1.7 


rur, 


= 55.8 , 


evC^/10-io = 2.7 



Table 2: Summary of x^(^^^) for the relevant data sets from the DAMA, CoGeNT and 
CRESST-II experiments. [Here Vrot = 200 km/s]. 

co-rotating halo, a 10% — 20% smaller effective Vrot value can easily be envisaged. There 
are also significant astrophysical uncertainties in this quantity. 
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Figure 2: DAMA, CoGeNT and CRESST-II favored regions of parameter space in the mirror 
dark matter model for Vrot = 200 km/s. 



The 2a {3a) favored region of parameter space is bounded by the contours where 
X^(^A',eVS^) = xfnin + 4(9). In figure 2 we plot the CRESST-II 2a favored region 
along with the 3a favored region of parameter space for DAMA and CoGeNT^. 



^Thc CoGeNT experiment does not discriminate electron recoils from nuclear recoils. In principle, 
e' inelastic scattering on bound electrons in the target can therefore also contribute to the signal at low 
energies Ej^ ~ keVee. Previous work[34, 11] has indicated that the e' contribution can be comparable 
to the rate of nuclear recoils. However, the previous analysis [34, 11] assumed that the flux of e' at the 
detector was much larger than the flux of mirror nuclei due to the large dispersion velocity vo[e'] ^ Vrot- 
In reality this cannot be the case. A larger flux of e' would lead to a greater rate of e' capture in the 
Earth c.f. mirror nuclei capture, and hence this would lead to an increasing mirror electric charge within 
the Earth, Q'^. In fact, we expect Q'^ to be generated such that the flux of e' and mirror nuclei at the 
Earth's surface approximately equalize. Taking this effect into account reduces the expected e' flux by 
more than an order of magnitude, leaving nuclear recoils as the dominant contribution to the rate. 



Experiment 



x'^{iTLin)/d.o.f. 



Best fit parameters 



DAMA (annual mod.) 5.5/10 

CoGeNT (spectrum) 15.8/12 

CRESST (spectrum) 0.3/3 



^rot 
Vrot 
Vrot 



210 km/s 
201 km/s 
250 km/s 



^/lO 
^/lO 



-10 = 3.1 

-10 = 1.6 
-10 ^ 1,7 



Table 3: Summary of x^(mm) for the relevant data sets from the DAMA, CoGeNT and 
CRESST-II experiments for A' = Fe'. 

Figure 2 indicates that in the high mass region, rriA' ~ 50mp including A' ~ Fe', 
DAMA encompasses lower ey/^ values. This is because dark matter scattering on Io- 
dine becomes kinematically favored and can dominate over Sodium scattering in this 
parameter region. The resulting lower values of ey/C,A' produce a favored region of pa- 
rameter space overlapping with that of the other experiments. If the spectrum of mirror 
dark matter elements has some resemblance to ordinary matter, we might expect Fe' to 
be the dominant heavy mirror metal component. To study this possibility further we 
fix A' = Fe' {rriA' = 55.8mp) and vary Vrot in the range: 160 < froi(km/s) < 250. The 
resulting Xmin values are given in table 3. Expanding around these Xmin values we obtain 
the CRESST-II 2a and the DAMA, CoGeNT 30" favored parameter regions in figure 3. 
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Figure 3: DAMA, CoGeNT and CRESST-II favored regions of parameter space in the mirror 
dark matter model for A' = Fe' (i.e. m^/ = 55.8mp). 



4 Two examples, PI and P2 

As discussed earlier, the systematic uncertainties are potentially quite significant for 
CoGeNT. The favored region can move up or down by ~ 20% as future data is collected 
and backgrounds are better understood. A small channeling fraction or quenching factors 
outside the range considered in Eq.(12), can also move the DAMA region somewhat. In 
this light, figures 2,3 offer some encouragement as they indicate some parameter space 
where all three experiments can be explained by Fe' scattering. At this point one could 
combine all three of these experiments and do a fit using the combined x^- However, 
in view of the potentially large systematic uncertainties it is not clear how useful such 
a combined analysis would be. Instead we consider two example reference points. The 
first one is located near the overlapping allowed regions indicated in figure 2: 



PI: A' = Fe {mpe' ^ 55.8mp), Vrot = 200 km/s, eJ^Pe' = 2.2 x 10 



-10 



(13) 



Figure 4a,b,c, compares the predicted rates assuming the PI parameter set with the 
relevant DAMA, CoGeNT and CRESST-II data. 
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Figure 4a: DAMA annual modulation spectrum for mirror dark matter with parameter choice 
PI (solid line). The separate contributions from dark matter scattering off Sodium (dashed 
line) and Iodine (dotted line) are shown. In this example q^a = 0.33, g/ = 0.20. 
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Figure 4b: CoGeNT spectrum for mirror dark matter with parameters PI (solid line). In this 
example qce = 0.17. 
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Figure 4c: CRESST-II spectrum for mirror dark matter with parameters PI (solid line). The 
signal component (dashed line) and background component (dotted line) are also shown. 
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For the example point PI, the CRESST-II signal is dominated by Fe' scattering on 
the Ca target element. Of course the reasons for this are purely kinematical being due 

to mpe' ~ rnca- 

The change in sign of the annual modulation amplitude at low energies indicated 
in figure 4a will not necessarily hold if there are lighter more abundant components 
such as O', Ne'. The positive contributions to the annual modulation amplitude from 
these lighter components can overwhelm the negative contribution due to the heavier 
components at low energies. We illustrate this point with the following example: 



P2:A' = Fe, eJ^Pe' = 1.7 x 10"^^ ^o'/^Fe' = 10, Vrot = 210 km/s . 



(14) 



Figure 5a,b,c, compares the predicted rates for the reference point P2 with the relevant 
DAMA, CoGeNT and CRESST-II data. 
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Figure 5a: DAMA annual modulation spectrum for mirror dark matter with parameters P2 
(solid line). The separate contributions from Fe' interactions (dashed line) and O' (dotted line) 
interactions are also shown. In this example qi\fa = 0.34, qj = 0.20. 
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Figure 5b: CoGeNT spectrum for mirror dark matter with parameters P2 (solid line). The 
separate contributions from Fe' interactions (dashed line) and O' (dotted line) interactions are 
also shown. In this example qce = 0.17. 
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Figure 5c: CRESST-II spectrum for mirror dark matter with parameters P2 (solid line). 
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Figures 5 show that with low Vrot ~ 220 km/s hght metal components (~ O') and 
lighter components (H', He') can be kinematically suppressed in DAMA, CoGeNT and 
CRESST-II. Future experiments at lower thresholds, such as Texono[35], should be able 
to probe such light components. 

Observe that the dark matter signal contribution in CRESST-II is relatively low in 
these examples. However it could easily be enhanced in many ways. Increasing e^/^i 
will increase the rate in CRESST-II and improve the fit in DAMA and is therefore an 
obvious possibility. Another possibility is that the CRESST-II energy scale might be 
overestimated by e.g. ~ 10%. If CRESST-II events are at lower energies then this can 
increase the expected rate in CRESST-II. Also, a slight increase in ffi of order ~ 20% 
can greatly improve the fit in CRESST-II. Such small changes of m do not significantly 
affect DAMA or CoGeNT since these experiments are sensitive to dark matter in the 
body of the Maxwellian distribution[ll] while CRESST-II is probing interactions in the 
tail of the dark matter velocity distribution (in this example). 
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Figure 6: DAMA spectrum for mirror dark matter with parameters PI (solid line), P2 (dashed 
line). In this example q^a = 0.34, qj = 0.20. 

Future data from DAMA, CoGeNT, CRESST-II and other experiments will obviously 
be able to constrain the parameter space within the mirror dark matter framework. As 
discussed recently [36], a particularly striking diurnal modulation signal should be ob- 
servable for a detector located in the southern hemisphere, and maybe even in a detector 
in the northern hemisphere at low latitudes, such as detectors in Jin-Ping Underground 
laboratory. In the meantime, we must rely on annual modulation and spectrum data. 
In figure 6 we give the predicted spectrum for DAMA/Libra and Figure 7 the predicted 
annual modulation spectrum in CoGeNT for each reference point, PI, P2 [detection 
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efficiency = 1 for these figures]. Clearly, the initial annual modulation amplitude mea- 
sured by CoGeNT to be A ^ 0.46 ± 0.17 cpd/kg/keVee averaged over the energy range: 
0.5 < i?(keVee) < 3.0, is much larger than that predicted by our example points. It 
will be interesting to see if CoGeNT's hint of a large annual modulation amplitude is 
conffimed by CDEX, C-4, and other experiments. 
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Figure 7: CoGeNT annual modulation spectrum for mirror dark matter with parameters PI 
(solid line), P2 (dashed line). In this example qce = 0.17. 

Let us make a few comments. Firstly, there are alternative explanations of DAMA, 
CoGeNT and CRESST-II within the mirror dark matter framework. For example, if 
Vrot is high (~ 240 km/s) then DAMA and CoGeNT can be mainly detecting the lighter 
metal components ~ O'-Si' [11, 12]. In this case, it is possible that CRESST-II might 
be seeing a less abundant heavier component (e.g. Fe'), or even the same component if 
Vrot is very high ~ 260 km/s. Secondly, more generic hidden sector dark matter can also 
explain the data in a similar way to the mirror dark matter case. One simply requires a 
hidden sector with two or more stable particles charged under an unbroken f/(l)' gauge 
interaction which is kinetically mixed with standard L'^(l)y[12, 10]. These ingredients 
can arise, for example, in spontaneously broken mirror models[37]. 

Finally, observe that the experiments constrain e^J^A' but do not allow for a precise 
determination of e or S,a' separately. It is difficult to estimate ^a'- Mirror metals might 
be generated in mirror stars, presumably at an early epoch[38]. It might also be possible 
to produce the metal component in the early Universe, although this is disfavored in 
the simplest scenarios with high reheating temperature [27]. A third possibility is that 
the average metal component is actually very low in the galaxy, {^a') ^ 10^^, but 
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its abundance is greatly enhanced in the galactic disk. One possible mechanism for 
this is the influence of the galactic magnetic field [39] which can potentially exclude (or 
suppress) light elements (m^' ~ ""^o') from entering the disk[24] (see also ref.[40]). This 
could thereby enhance the proportion of heavier components whose abundance should 
increase until the pressure equalizes with that of the rest of the halo (at the same radial 
distance). In this case C,a' ~ 1 is possible within the disk, i.e. at the detector's location. 
The above possibilities suggest a broad range for ^a''- 10~^ ~ C,a' ~ 1- If the direct 
detection experiments are due to mirror dark matter with e^/^ ~ few 10~^° then this 
suggests that e is in the range^" ~ 10~^° — 10~*. The most stringent laboratory limit on 
e arises from invisible decays of orthopositronium[43, 44], e ~ 1.5 x 10~^. An important 
proposal exists [45] for a more sensitive orthopositronium experiment which can cover 
much of the e range of interest. Also, note that values of e ~ 10^^ can also be probed 
in early Universe cosmology [46]. In particular, forthcoming results expected from the 
Planck mission might shed some light on e in the near future. There are also other 
interesting possible effects of the kinetic mixing[15, 47]. 

5 Constraints from CDMS and XENONIOO 

Let us now examine the question of the compatibility of this model with the constraints 
from the XENON100[48], CDMS/Ge[49] and CDMS/Si[50] experiments^^ This ques- 
tion is somewhat non-trivial since the answer depends sensitively on the recoil energy 
threshold which (typically) has at least 20% uncertainty and often the subject of con- 
troversy. The XENONIOO experiment, which has the largest exposure, has especially 
poor knowledge of energy calibration in the important low energy region (see ref. [54, 55] 
for relevant discussions). In view of the above, we explore the compatibility question by 
estimating the energy threshold for which the parameter point PI can be excluded at 
95% C.L. [The constraint on P2 is similar to that of PI]. We have taken into account 
the relevant detection efficiencies, exposure time etc, although we conservatively ignore 
effects of detector resolution for the XENONIOO experiment. Our results are given in 
table 4. Note that the 95% C.L. limits in the second column arise since CDMS/Ge, 
CDMS/Si and XENONIOO observed 2, 0, 1 events respectively in the low energy region. 
Table 4 indicates that mirror dark matter with parameters given by the example PI 
is consistent with the CDMS/Ge, CDMS/Si and XENONIOO data provided that the 
energy scale has been underestimated by of order 20 — 30% respectively. This level of 
energy scale uncertainty is within critical assessments of these experiments discussed 
in the literature. We conclude that the considered mirror dark matter explanation of 
the DAMA, CoGeNT and CRESST-II data is consistent with the data from CDMS/Si, 
CDMS/Ge and XENONIOO experiments when reasonable systematic uncertainties in 



^^Siniilar values of e can be motivated from an astrophysical argument. Mirror dark matter is dissi- 
pative and radiative cooling would cause the halo to collapse on a time scale of a few hundred million 
years unless a substantial heat source exists. Kinetic mixing induced processes in the core of ordinary 
supernova (such as e'e' production via plasmon decay[41]) can supply the energy needed to compensate 
for the energy lost from the halo due to radiative cooling if e ~ 10~^[23, 42]. 

^^There are also lower threshold analysis by the XENON10[51] and CDMS collaborations[52]. However 
it has been argued [53] that neither analysis can exclude light dark matter when systematic uncertainties 
are properly taken into account. 
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Experiment 


95% C.L limit 


Eg^* threshold 


^nom threshold 


100 


fpnom 




CDMS/Ge 

CDMS/Si 

XENONIOO 


6.3 
3.0 

4.7 


13.1 keV 
8.5 keV 
12.9 keV 


10.0 keV 
7.0 keV 
8.4 keV 


23.7% 
17.6 % 
34.9 % 



Table 4: Estimated energy threshold (Ef *) for which the CDMS/Ge, CDMS/Si and 
XENONIOO experiments are consistent with the mirror dark matter expectations assum- 
ing the example point, PI. E^"^ is the nominal energy threshold of each experiment. 
The last column gives the percentage at which the energy threshold needs to have been 
underestimated for the point PI to be consistent with these experiments at 95% C.L. 

energy scale are included. It should also be noted that the two events seen in CDMS/Ge 
and also the events seen just above the threshold in the Edelweiss experiment [56] are 
both compatible with dark matter interactions within this model given these energy scale 
uncertainties. 

Finally, the KIMS experimental] is currently running with a ~ 100 kg Csl target. 
They have reported a (2o") limit on the rate of dark matter interactions from a Pulse 
Shape Analysis (PSA) of their spectrum: R < 0.02 cpd/kg/keVee for 3.0 < EjiikeVee) < 
4.0. We find that the Iodine interaction rate in a Csl target predicted by mirror dark 
matter for the particular point PI is i? ~ 0.006 cpd/kg/keVee in the 3 < EnikeVee) < 4 
bin and falling sharply at higher energies. [Assuming here that KIMS has the same 
resolution and qj quenching factor as the DAMA/Libra experiment]. Thus, their PSA 
does not constrain mirror dark matter. However, the KIMS experiment also aims to 
measure/constrain the annual modulation amplitude. The annual modulation results 
are anticipated to be reported within the coming year. Assuming that they can reach 
a similar energy threshold to DAMA, they can potentially see a positive signal for a 
signficant portion of parameter space (c.f dotted line in figure 4a). 

6 Conclusion 

In conclusion, we have examined the new results from the CRESST-II experiment along 
with the latest results from DAMA and CoGeNT experiments in the context of the mirror 
dark matter framework. In this framework dark matter consists of a spectrum of mirror 
elements: H', He', O', Fe', ... of known masses. Under the simplifying but reasonable 
assumption that DAMA, CoGeNT and CRESST-II might be observing a particular dark 
matter component. A', we have found that mirror dark matter models can explain the 
data from each experiment. In particular, we have found that each experiment can 
be explained by A' ~ Fe' interactions if e^/^pe' ~ 2 x 10^^" and Vrot ~ 200 km/s. 
Other regions of parameter space are possible. We have also shown that the considered 
explanation is consistent with the results of the other experiments when reasonable 
systematic uncertainties in energy scale are considered. 
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